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The electronic structure and magnetism of chromia (corundum- type Cr2 03) are studied using full- 
potential first-principles calculations. The electronic correlations are included within the LSDA-I-U 
method. The energies of different magnetic configurations are very well fitted by the Heisenberg 
Hamiltonian with strong exchange interaction with two nearest neighbors and additional weak in- 
teraction up to the fifth neighbor shell. These energies are insensitive to the position of the oxygen 
states, indicating that magnetism in Cr203 is dominated by direct exchange. The Neel temperature 
is calculated using the pair-cluster approximation for localized quantum spins of magnitude 3/2. 
Very good agreement with experiment is found for all properties including the equilibrium volume, 
spectral density, local magnetic moment, band gap, and the Neel temperature for the values of U 
and J that are close to those obtained within the constrained occupation method. The band gap is 
of the Mott-Hubbard type. 



Corundum-type Cr203 is one of the antiferromagnetic 
transition-metal oxides which present a significant chal- 
lenge for electronic band theory due to the correlated 
character of the partially filled, spin-polarized 3d shell. It 
also has numerous applications in electronic devices, fuel 
cell electrodes, gas sensors, heterogeneous catalysis, and 
thermal barrier coatings. Surface properties of Cr203 are 
of particular interest. The nominally polar, but compen- 
sated (0001) surface exhibits structural phase transitions 
which are poorly understood fii^ as well as an uncompen- 
sated surface moment'^ which may be useful in spintronic 
applications. It is therefore very desirable to establish 
whether electronic correlations can be reliably included 
in first-principles calculations in a way that would accu- 
rately predict structural, electronic, and magnetic prop- 
erties. 

As expected for a transition-metal oxide, conven- 
tional density-functional theory (DFT) studies of bulk 
Cr203'*'^'^ have shown that local spin density approxima- 
tion (LSDA) or the generalized-gradient approximation 
(GGA) for the exchange-correlation potential are unable 
to reproduce the electronic and magnetic properties of 
bulk Cr203. Rohrbach et air- performed a GGA-fU cal- 
culation for Cr203 using the simplified (spherically av- 
eraged) U — J correction^ and obtained more reasonable 
results for the band structure. However, this approach is 
inaccurate for structural and magnetic properties. First, 
as is typical for transision-metal compounds, both GGA 
and the LSDA+U corrections reduce the LSDA overbind- 
ing problem in Cr203 and increase the equilibrium vol- 
ume which becomes close to the experimental value. 
However, the use of both GGA and Hubbard U in the 
GGA-I-U method^ results in an overcorrection, so that 
the equilibrium volume becomes 7% too large. Also, 
the magnetic energies found by Rohrbach et al. within 
this method are incompatible with the experimental Neel 
temperature of about 308 K; they are too small roughly 
by a factor of 5. Mosey et al.^ calculated the U and J 
parameters using the unrestricted Hartree-Fock method 
and found U — J = 7.7 eV, which is quite obviously too 



large for Cr, as is typical for Hartree-Fock-like methods. 
They also studied the structural and electronic proper- 
ties of Cr2 03 using the spherically averaged LSDA-I-U 
implementation* and found that the structural proper- 
ties and the band gap come out best at U — J — 4 eV. 
For this value of U — J the magnetic moment is some- 
what reduced from the "ideal" (ionic) value of 3 fiB due 
to hybridization with oxygen, which is also in agreement 
with experimental dataJ^ 

In this paper we analyze the properties of Cr203 in 
more detail with a particular emphasis on its magnetism. 
We found that the LSDA-I-U method can provide very 
good agreement with experiment for structural, spectral, 
and magnetic properties with the same values of U and 
J which are close to those found from the constrained 
occupation method within DFT. 

The first-principles calculations were carried out using 
the projected augmented wave (PAW) methodiiii^ im- 
plemented within the Vienna ab initio simulation package 
(VASP)ji^'i^ The 2s states on O were treated as valence 
states. We used the rhombohedral primitive cell for the 
corundum structure in all calculations except those in- 
volving complicated magnetic configurations (see below). 
The Monkhorst-Pack scheme^^ based on the 4 x 4 x 4 fc- 
point grid was employed for the Brillouin zone integrals, 
which were calculated using the tetrahedron method with 
Blochl corrections. The plane- wave energy cutoff was 520 
eV. These parameters ensured the total energy conver- 
gence to 2 meV/atom. Densities of states (DOS) were 
calculated using the 8x8x8 fc-point grid. 

We employ the LSDA+U method in its full spherically 
symmetric formJ^ This extension is important for Cr203 
where correct representation of crystal field and exchange 
splittings within the partially filled 3d shell is critical. 
Surface energetics are even more sensitive to the cor- 
rect treatment of these unfilled shells; in fact, we found 
that the errors introduced by the spherically averaged 
LSDA-I-U ansatz are intolerable for the Cr2 03(0001) 
surfaceJ^ The double-counting term is taken in the fully 
localized limit J^ii^ 
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Reasonable values of U and J can usually be obtained 
within DFT using the constrained occupation method 
These values are better suited for use in LSDA+U cal- 
culations compared to the Hartree-Fock values, because 
they include the screening of the id shell by the remaining 
electrons. We calculated U and J in this way using the 
full-potential linear augmented plane-wave (FLAPW) 
method implemented in the FLEUR package^SS. In these 
calculations the GGA approximation was used. We took 
the rhombohedral primitive cell of Cr203 containing four 
Cr atoms and set all the structural parameters to their 
experimental values 3i The 3d electrons on one or two 
(most distant) Cr atoms were formally treated as open 
core shells, i.e. an integer occupation of these orbitals (for 
each spin projection) was enforced, and their hybridiza- 
tion with all other electrons was turned off. The U and 
J parameters are then found by comparing the LSDA 
total energies for different charge and spin occupations 
of the id orbital(s) with their Hartree-Fock expressions 
(the latter are equal to the "double-counting" terms in 
LSDA-hU). 

The constrained occupation method is somewhat am- 
biguous because the U parameter depends on the charge 
state of the ionS^. (On the contrary, the J parameter is 
usually very well defined; we found this to be the case for 
Cr203 as well.) Although the formal charge state of the 
chromium ion is Cr'^"'", we find U and J with respect to 
the id^ state. The reason is that the screening properties 
of the valence electrons depend primarily on the charge 
density distribution in the crystal, which is typically very 
close to the superposition of atomic charge densities. In- 
deed, the formal occupancy of the Cr id orbital within 
the 2.5 a.u. muffin tin sphere is about 4.2 in FLAPW. 
Specifically, the exchange parameter J — 0.58 eV was 
found by treating 3d electrons on one Cr atom as core and 
considering the energy difference between the id^id^^^, 

3(i|3dJ and id'^id^^^ configurations. The Hubbard param- 
eter U = 3.3 eV was found by treating the id shells on the 
two most distant Cr atoms as open cores and consider- 
ing the energy difference between the id\id2 and id\id2 
configurations, where 1 and 2 refer to the two different 
sites (the contribution of J to the energy differences was 
subtracted). Since the total number of electrons in the 
cell is the same for both configurations, there is no need 
to include the Fermi level correction. 

In the following calculations we fixed J at its calcu- 
lated value of 0.58 eV and varied U . The electronic struc- 
ture is calculated in the ground antiferromagnetic state 
with relaxed structural parameters (commonly denoted 

as H 1 — in accordance with the ordering of Cr spins 

along the z axis of the rhombohedral cell). The depen- 
dence of the equilibrium volume, magnetic moment, and 
band gap on U is shown in Fig.[TJ We see that the volume 
and the band gap agree quite well with experimental data 
at [/ = 4 eV. Other structural parameters also agree with 
experiment. The calculated angle between the rhombo- 
hedral unit vectors is 55.11° compared to the measuredSi 
angle of 55.13°. The shortest distance between Cr atoms 
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FIG. 1: Atomic volume (a), magnetic moment (b), band gap 
(c) and Neel temperature (d) as functions of the Coulomb 

U parameter for the antiferromagnetic H h— state. The 

value of J is fixed at 0.58 eV. The horizontal lines denote 
exper imental values; those in panel (b) are from both Refs. 



along the [111] axis is 2.646 A vs the measured 2.650 A. 

The magnetic moment aXU — 4.0 eV is 2.86 /is, i.e. it 
is somewhat reduced compared to the "ideal" ionic value 
of 3 /is corresponding to a fully localized spin 3/2. Ex- 
perimentally, the most recent neutron polarimetry mea- 
surement gives the sublattice magnetization of 2.48 ij^et^, 
which is notable lower compared to an older estimate of 
2.76 ^B-^'^ The magnetic moment is smaller than 3 /is 
due to two effects: (1) hybridization with oxygen, which 
is included in our calculation, and (2) the quantum "zero- 
point spin deviation," which is absent in DFT. The zero- 
point deviation in Cr2 03 was estimated^^ to be about 
8%, which amounts to 0.24 ^b- Thus, keeping in mind 
the uncertainties related to the definition^ measurement, 
and calculation of the magnetic moment, its calculated 
value at [/ w 4.0 eV is completely reasonable. Of course, 
the poorly defined reduction from i jiB can not be used 
as an indicator of the quality of agreement with exper- 
iment. We also note that the local magnetic moment 
depends very weakly on the magnetic configuration of 
Cr203. In the ferromagnetic state the local moment 
within the muffin-tin sphere is 2.94 /i^, while the magne- 
tization is exactly equal to 3 /is per Cr site, as expected. 
The orbital moment in the calculation with spin-orbit 
coupling is small (about 0.04 hb) and antiparallel to the 
spin moment, in agreement with the experimental^ g- 
factor of 1.97 and with the general rule for atomic shells 
that are less than half filled. 
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FIG. 2: Orbital projected DOS for the antiferromagnetic 

H 1 — ground state of bulk Cr2 03, calculated within the 

LSDA+U with U = 4.0 eV, J = 0.58 eV. (a) 0-2p, (b) Cr-eg 
and (c) Cr-t2g. Majority and minority-spin DOS are plotted 
with different signs. 



The band gap of 3.07 eV is somewhat smaller than the 
experimental value of 3.4 eV, but greater than that found 
in Ref. 't' using GGA+U with U - J = A eV. Underesti- 
mation of the addition energy is a common feature of the 
LSDA-I-U method, which is well known, for example, for 
Gd and other 4/ elements. Further, the density of states 
shown in Fig. [2] is in excellent agreement with X-ray pho- 
toemission data^ii^. Namely, the sharp and narrow peak 
at low binding energies separated by a (pseudo)gap from 
the rest of the valence band is very well reproduced. Fig. 
[5] also shows the partial DOS decomposition into O and 
Cr contributions from t2g and eg states. (The t2g and 
eg subbands are well defined because the ligand field of 
the Cr site is approximately octahedral.) As seen from 
Fig. [21 the peak at low binding energies corresponds to 
the filled Cr t2g spin subband with some admixture of 
oxygen p orbitals. At least some of this admixture is fic- 
titious, because the Cr d orbitals extend into oxygen's 
projection spheres. (The same ambiguity is involved in 
the definition of the Cr magnetic moment.) Since non- 
magnetic Cr2 03 is metallic with the Fermi level lying 
inside the t2g subband which is separated by a gap from 
the oxygen p band, the insulating gap forms by the Mott- 
Hubbard mechanism. 

Thus, the available structural and spectral properties 
of Cr203 are well reproduced by the spherically symmet- 
ric LSDA+U method with U « 4.0 eV and J = 0.58 eV. 
This value of U is somewhat larger than that given by the 
constrained occupation method, but the latter was found 
using FLAPW calculations with a different muffin-tin ra- 
dius. 

We now focus on the magnetic energetics which 
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TABLE I: The exchange parameters J„ fitted using the total 
energies of 12 magnetic configurations calculated for different 
values of U with J fixed at 0.58 eV. Each pair of sites con- 
tributes JijGiej to the total energy, where is the unit vector 
parallel to the direction of the local moment at site i. The 
last column A shows the mean-square misfit in the fitting of 
the total energy. J„ and A are given in meV. 
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provides another stringent test of the validity of the 
LSDA-fU method for Cr203. We calculated the total 
energies of 12 different magnetic configurations includ- 
ing the ferromagnetic state, three simple antiferromag- 
netic orderings (-1 1 — , + -\ , and 4 h), and 

eight additional, arbitrarily chosen spin configurations in 
the hexagonal unit cell of the corundum structure which 
includes six formula units. The ground antiferromagnetic 

state has H 1 — ordering in agreement with experiment; 

its relaxed structure was fixed for other magnetic order- 
ings. The calculated total energies are well fitted by the 
conventional Heisenberg Hamiltonian with exchange in- 
teraction in five coordination spheres. The exchange pa- 
rameters fitted by the least squares method for different 
values of U are listed in Table [J The first, second and 
fifth nearest-neighbor interactions are antiferromagnetic, 
whereas the third and fourth ones are ferromagnetic. All 
exchange parameters except the third one reinforce the 
ground state. Also, we see that the first and second ex- 
change parameters are significantly larger than those for 
more distant neighbors; this behavior is natural for an 
insulator. In Cr203 the parameter Ji corresponds to the 
short bond along the z axis, and J2 to the nearest neigh- 
bor bond within the buckled Cr layer. Each Cr atom has 
one Ji neighbor and three J2 neighbors. As the value of 
U is increased, all antiferromagnetic exchange parame- 
ters systematically decrease, while the ferromagnetic pa- 
rameters systematically increase. 

Let us compare our results with those derived from 
the spin wave dispersion measured by inelastic neutron 
scatteringi^ The exchange parameters fitted to those 
dispersion curves are also dominated by antiferromag- 
netic Ji and J2. However, more distant neighbors in 
that fitting^^ are much less important. For example, 
J5/J1 « 1/40 in Ref. HI 1/7 found here. Also, there 
is a notable disagreement in the J2/J1 ratio: Ref. 
found J2/J1 ~ 0.45 versus 0.76 for [/ = 4 eV in Table [J 
We have verified that the anisotropy J2/J1 is almost un- 
changed if only the lowest-lying spin configurations are 
included in the fitting. 

We also considered the possibility that the exchange 
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parameters may be affected by lattice distortion below 
the Neel temperature due to magnetostructural coupling. 
In order to study this possibility, we need to know how 
the exchange parameters depend on the structural pa- 
rameters, and how the latter change between room tem- 
perature and liquid nitrogen temperature where the spin 
wave spectrum was measured. We found that the values 
of Ji and J2 are very sensitive to the lattice parameters; 
a 1% increase or decrease in the first or second neighbor 
bond length leads to a 25-50% decrease or increase of the 
corresponding exchange parameter (this was roughly es- 
tablished by varying the Bravais lattice parameters of the 
rhombohedral cell and fitting the four simple magnetic 
orderings to the Heisenberg model). Thus, the J2/J1 ra- 
tio is very sensitive to the c/a ratio (or any lattice distor- 
tion that changes the ratio of the first and second bond 
lengths R2/R1). 

To our knowledge, no experimental data are available 
on the thermal contraction of Cr203 below room temper- 
ature. In order to assess the degree of magnetostructural 
coupling, we have fully relaxed the structure for all the 12 
magnetic configurations that were included in the fitting 
of exchange parameters. The magnetostructural coupling 
is, in general, quite appreciable; the R2/R1 ratio varies 
between the minimum of 1.072 in the H h config- 
uration (where the nearest-neighbor pairs are all paral- 
lel, and next-nearest antiparallel) and the maximum of 

1.103 in the + H configuration (where the situation 

is reversed). The overall trend, as expected in general 
for antiferromagnetic coupling, is for each of these bonds 
to shorten when the corresponding spins are antiparallel 
and lengthen when they are parallel. Although this ef- 
fect seems to be rather large, the actual H 1 — state 

has both first and second-neighbor pairs antiparallel, so 
that both should lengthen somewhat in the paramagnetic 
state. Among our 12 configurations, 5 have the same ra- 
tio Np'^ /N^p of the number of parallel and antiparallel 
pairs in the i-th coordination sphere for i = 1 and « = 2, 
i.e. these 5 configurations approximately represent the 
change of structure as a function of the antiferromagnetic 
order parameter. The R2/R1 ratio for these 5 structures 

varies between 1.088 in the H 1 — state and 1.084 in 

the ferromagnetic state. The 3 intermediate states, all 
of which have twice as many parallel pairs than antipar- 
allel ones, all have R2/R1 — 1.087. The bond lengths 

themselves change by about 0.7% between the H \ — 

state and these three intermediate states. This variation 
is obviously too small to explain the disagreement with 
experiment in the J2/J1 ratio. 

Thus, in spite of the overall agreement with experiment 
for most properties, a discrepancy in the J2/J1 ratio re- 
mains. It is possible that our Heisenberg model fitting is 
not fully applicable to small deviations from the ground 
state, for which the linear response technique should be 
a better fit. Further, spin waves are more sensitive to 
distant couplings compared to the energy fits or thermo- 
dynamic properties, which makes the fits from spin wave 
data and from the overall energetics statistically inequiv- 



alent. Finally, many-body effects beyond LSDA-I-U may 
play a role in renormalizing the exchange parameters. 
We also note that impurities and thermal spin disorder 
should be more effective in destroying distant couplings 
compared to nearest-neighbor ones; this may explain the 
larger role of couplings beyond J2 in our fitting compared 
to experimental spin wave results. 

The mechanisms of exchange interaction in Cr203 are 
not well understood. Complicated crystal structure with 
the presence of many electronic orbitals of different sym- 
metry and many Cr-O-Cr links at different angles greatly 
complicates the empirical analysis. Both superexchange 
and direct exchange interactions have been invoked to 
explain the magnetic structure of Cr203i^i2£i^iH Direct 
exchange interaction may be interpreted in terms of hop- 
ping of electrons from one Cr ion to another across the 
insulating gap; the energy denominator involved in this 
process is the Mott-Hubbard splitting. The hopping can 
only be effective between the orbitals that are able to hy- 
bridize. From the DOS decomposition shown in Fig. [2] it 
is clear that the t2g orbitals very weakly hybridize with 
Eg orbitals on the neighboring Cr atoms. Therefore, the 
contribution of Cg orbitals to direct exchange can be ne- 
glected. The t2g subband is split off by crystal field and 
exactly half-filled, therefore direct exchange should be 
antiferromagnetic. The superexchange involves hopping 
between Cr and O ions; the energy denominator involves 
the gap between the oxygen p states and the unoccupied 
Cr states. 

In order to reveal the mechanism responsible for ex- 
change interaction in Cr203, we use the following trick. 
A fictitious external potential V is coupled to the oxygen 
p orbitals, which adds the term Ey = ^TrnJ^^j, to the 
total energy, where the trace is taken over orbital and 
spin indices, and n'^^, is the density matrix of the oxy- 
gen p states defined inside the muffin tin sphere of 1.2 
a.u. This density matrix is calculated using the standard 
LSDA-I-U machinery. When V is set to a negative value, 
the filled oxygen p states are pushed down to lower ener- 
gies, which suppresses superexchange, but not direct ex- 
change. Weak hybridization between the filled t2g states 
and the oxygen p states pushes them apart at = 0. 
When the p states are pushed down by increasing V, this 
repulsion is reduced and the t2g states also move down, 
thereby increasing the band gap. Since direct exchange is 
sensitive to the band gap, for better comparison we com- 
pensate this increased band gap by reducing U on Cr. 
This is done in such a way that the distance between the 
center of mass of the filled t2g band and the conduction 
band minimum (CBM) remains the same as at F = 0. 

Starting from the state with U = 3.5 eV, we added 
V of —12 eV and —24 eV and calculated the energies 
of four magnetic states in the rhombohedral primitive 
cell of Cr203. These calculations were performed using 
the FLAPW methodi^ the results are listed in Table 
im (Those for y = are within 10% of VASP results.) 
The self-consistent downward shift of the p states is much 
smaller than Vrimm because V is strongly screened. The 



5 



TABLE II: Energies of three simple magnetic configurations 

relative to the ground H h— state (in meV per formula unit) 

as a function of the fictitious external potential V applied to 
the oxygen p orbitals (in eV). Reduced values of U on Cr 
compensate for the increased band gap (see text for details). 
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distance A from the oxygen p band maximum to CBM 
is also listed in the table. 

The results listed in Table |TT] clearly show that the mag- 
netic energies are insensitive to the position of the oxygen 
p band. The reduction of magnetic energies produced by 
adding negative V at constant U is due to the fact that 
the t2g band gap increases due to dehybridization from 
oxygen. Once U is decreased to bring the band gap to 
its original value, the magnetic energies are essentially 
unchanged compared to their values at y = 0; in fact, 
they even increase somewhat. On the other hand, we've 
seen above that the magnetic energies are very sensitive 
to the value of U which is responsible for the band gap. 
This behavior leads us to a striking conclusion that, con- 
trary to the common belief, superexchange plays no role 
in magnetism of Cr203. Antiferromagnetism is due ex- 
clusively to direct exchange, which, as mentioned above, 
is antiferromagnetic because the magnetically active t2g 
subband is half-filled. It is likely that superexchange in 
Cr203 is highly ineffective because the Cr-O-Cr angles 
are close to 90°, while the overlap between O states and 
Cr t2g states is small. On the other hand, the overlaps 
between t2g states on neighboring Cr atoms are quite 
large; the t2g bandwidth in the ferromagnetic state at 
V = —24 eV is approximately 1.5 eV. 

We now calculate the Neel temperature T^. We saw 
above that the local moments on Cr atoms are well local- 
ized, and the energies of spin configurations are well rep- 
resented by the Heisenberg spin Hamiltonian. We there- 
fore adopt the quantum Heisenberg model for localized 
spins of magnitude 3/2 with the exchange parameters 
listed in Table HI Since each spin is strongly coupled 
only to four neighbors (one with Ji and three with J2), 
the mean-field approximation can not be reliably used. 
However, the antiferromagnetic interaction is not frus- 
trated. The important spin correlations should be gen- 
erated by the dominant exchange interaction with four 
nearest neighbors. The network of bonds corresponding 
to Ji and J2 is very weakly connected; the shortest closed 
path on this network is 6 bonds long. Therefore, it is suf- 
ficient to capture the pairwise spin correlations. In this 
situation, the pair cluster approximation appears to be 



an obvious choice. This approach is a special case of the 
cluster variation method when the set of clusters includes 
only pairs of sites. Here we follow the formulation of Ref. 

which can be directly applied to our case. The details 
of this technique are described in the Appendix. 

The calculated is shown in Fig. [1] as a function of 
U. We see that the best agreement with experiment for 
T/v is obtained at the same value of J7 « 4 eV as for the 
structural and spectral properties explored above. This 
overall agreement is a strong indication that the essential 
details of the electronic structure of Cr203 are very well 
captured by the LSDA+U approximation. Physically, 
this success of LSDA+U is explained by the presence 
of fully filled and empty subbands which are strongly 
split by crystal and exchange fields; LSDA-fU usually 
reproduces such closed atomic-like subshells very well. 

In conclusion, we found that the spherically symmet- 
ric LSDA-fU method provides a very good description of 
structural, spectral, and magnetic properties of chromia 
with J = 0.58 eV found from the constrained occupation 
method and J7 « 4 eV, which is also close to the calcu- 
lated value. We found that the magnetic energies are well 
represented by a Heisenberg model with strong exchange 
interaction with nearest neighbors both in the plane and 
along the z axis and much weaker interaction with more 
distant neighbors. The artificial downward shift of the 
filled oxygen p states has almost no effect on the mag- 
netic energies, which proves that direct exchange is the 
dominant mechanism of magnetic interaction. 

This work was supported by NRC/NRI supplement to 
NSF-MRSEC and by the Nebraska Research Initiative. 
K. B. is a Cottrell Scholar of Research Corporation. 



APPENDIX 

Here we describe the application of the pair cluster 
approximation to Cr203 following a similar formalism of 
Ref. [33!. The energy of a quantum Heisenberg magnet 
(per unit cell) can be written as 

E=-IY1 ^^n^3J^3{^^ ' Sj) ^ ^^B,{St) , (A.l) 

where the summations are over the inequivalent sites in 
the unit cell. Si are quantum spin operators, rrii is the 
number of sites of type i in the unit cell, the num- 
ber of neighbors of site i that are of type j, and Bi the 
external magnetic field applied to site type i. All the 
thermodynamic properties can be obtained from the free 
energy which may be calculated by integrating the Gibbs- 
Helmholtz relation: 

f E{l3')df3', (A.2) 

To proceed we need to find the expectation values ap- 
pearing m Eq. dXTJ. In the pair-cluster approximation 
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they are calculated by introducing one- and two-site clus- 
ters with the following cluster Hamiltonians: 



HI 



TTl] _ 

-"2 — 



JijSi ■ Sj 



h^^Sf-h}:>S 



(A.3) 



where hi = Bi + (pi is the one-site "cluster field," h. 



(i) 



is the cluster field at site i for the pair cluster 



(ij). The one-site and two-site cluster fields are related 

through 4>i = J2j ''^ij'I^^P ■ The quantities (fif^ are treated 
as variational parameters and found by minimizing the 



free energy, i.e. requesting that dF/dcj), 



U) 



0. It can 



be shown that this variational condition ensures that the 
expectation value {Sf) is the same in all one-site and 
all two-site clusters containing site i. The expectation 
values (Si ■ S^) are calculated from the pair cluster (ij). 
Performing the integration in (|A.2p . we find 



F — — raiUij In Z!^ + — (n^ — 1) In , 



213^ 



1 



(A.4) 

where rii — riij is the total number of neighbors of site 

i, while Z\ = Trexp(-/3iJi) and Z'^ = Tr exp(-/3i/*^) 
are the one-site and two-site cluster partition sums. Eval- 
uation of Z\ is trivial; to find Z^ one first needs to di- 
agonalize H2 ■ 



Our goal here is to find the transition temperature; 

(7) 

therefore we may assume that all hf are small. The free 
energy is developed in these parameters, and then the 
parameters 0^"'^ are found by requiring that the variation 
of the free energy F vanishes. (The resulting equations 
are too cumbersome to be included here.) 

For Cr2 03 we assume the actual magnetic ordering 

H 1 — . All Cr sites are related by magnetic symmetry, 

which reduces the number of independent variational pa- 
rameters. The transition temperature is found by setting 
S = 3/2 and searching for the pole of the magnetic sus- 
ceptibility which is found from the one-site cluster: 



_9(5,f) _ 5(5 + 1) 



The resulting equation for the transition temperature 
has two solutions. The greater solution is Tpf, while the 
lower one is the fictitious "anti-Neel point" below which 
(5f ) = 0. The existence of an anti-Neel point for an- 
tiferromagnets is a well-known shortcoming of the pair- 
cluster approximation ;'^^'^^ which fails at low tempera- 
tures. However, in our case the anti-Neel point is much 
smaller than T/v which indicates that the pair cluster ap- 
proximation has a wide range of validity. Therefore, this 
method is expected to provide a very good approximation 
for Tat. 
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